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Abstract. This paper presents the theoretical basis of the 
fireball/blast wave model, and some implications of recent 
results on GRB source models and cosmic-ray production 
from GRBs. BATSE observations of the prompt 7-ray lumi- 
nous phase, and Beppo-SAX and long wavelength afterglow 
observations of GRBs are briefly summarized. Derivation 
of spectral and temporal indices of an adiabatic blast wave 
decelerating in a uniform surrounding medium in the limit- 
ing case of a nonrelativistic reverse shock, both for spheri- 
cal and collimated outflows, is presented as an example of 
the general theory. External shock model fits for the after- 
glow lead to the conclusion that GRB outflows are jetted. 
The external shock model also explains the temporal dura- 
tion distribution and clustering of peak energies in prompt 
spectra of long-duration GRBs, from which the redshift de- 
pendence of the GRB source rate density can be derived. 
Source models are reviewed in light of the constant energy 
reservoir result of Frail et al. that implies a total GRB en- 
ergy of a few x 10 51 ergs and an average beaming fraction of 
rs 1/500 of full sky. Paczyriski's isotropic hypernova model 
is ruled out. The Vietri-Stella model two-step collapse pro- 
cess is preferred over a hypernova/collapsar model in view 
of the X-ray observations of GRBs and the constant energy 
reservoir result. Second-order processes in GRB blast waves 
can accelerate particles to ultra-high energies. GRBs may be 
the sources of UHECRs and cosmic rays with energies above 
the knee of the cosmic ray spectrum. High-energy neutrino 
and 7-ray observations with GLAST and ground-based 7-ray 
telescopes will be crucial to test GRB source models. 



1 Introduction 

In only ten years, the study of gamma-ray bursts has grown 
from a narrow specialty involving a small group of high- 
energy astronomers to a mainstream field of research that 
impacts nearly all topics in astronomy. This change resulted 
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primarily from observational advances made with the Burst 
and Transient Source Experiment (BATSE) on the Compton 
Gamma Ray Observatory, and the discovery of X-ray after- 
glows and long wavelength counterparts made possible with 
the Beppo-SAX experiment. The BATSE results showed that 
the angular distribution of GRBs on the sky is isotropic, and 
that the GRB size distribution exhibits a strong flattening for 
faint bursts (Meegan et al., 1992). This behavior follows 
from a cosmological origin of GRB sources, with the de- 
cline in the number of faint bursts due to cosmic expansion. 
Follow-up X-ray observations with the Narrow Field Instru- 
ments on Beppo-SAX has permitted redshift determinations 
that firmly establish the distance scale to the sources of ^ 2 
s duration GRBs, which are those to which Beppo-SAX is 
sensitive. 

The redshifts of nearly 20 GRBs are now known, with the 
mean redshift near unity and the largest measured redshift 
at z = 4.5. The corresponding distances imply apparent 
isotropic 7-ray energy releases in the range from sa 10 51 - 
10 54 ergs. Recent results suggest that GRB emissions are 
strongly beamed, so that the total energy release is actually in 
the neighborhood of 10 51 ergs, corresponding to the typical 
total energies released in the kinetic outflow of a supernova. 
Delayed reddened enhancements detected in the optical light 
curves of a few GRBs could also be a consequence of a su- 
pernova emission component. If the beaming results are cor- 
rect, then many more sources of GRBs exist than are implied 
through direct statistical studies of detected GRBs. The im- 
plied rate of both aligned and misdirected GRB sources be- 
gins to approach the expected rate of Type Ib/c supernovae 
(SNe). These and other lines of evidence indicate that GRBs 
are related to a type of supernova. 

This review of the theory of GRBs briefly touches on some 
of the central questions in contemporary GRB studies. Af- 
ter summarizing the prompt 7-ray luminous and delayed af- 
terglow observations of GRBs, detections of X-ray emission 
line and absorption features, and some related phenomenol- 
ogy, the redshift distribution is presented. These redshifts are 
obtained from absorption lines in the optical afterglow, spec- 
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troscopy of the host galaxies and, in a few cases, from X-ray 
line observations. Given the redshift, the directed energy re- 
leases per unit solid angle can be derived. The theoretical 
basis for the fireball/blast wave model is outlined, and nu- 
merical calculations derived from the model are presented. 
According to the blast-wave model, a central source ejects 
an expanding relativistic shell of particles. Accelerated elec- 
trons radiate synchrotron photons to produce the long wave- 
length X-ray, optical, and radio afterglows, as well as the 
soft 7-ray emission emitted in the prompt phase (although 
arguments have also been made for a Compton origin of this 
emission). By contrast, a synchrotron self-Compton origin 
may explain the > 100 MeV emission detected with EGRET 
from a handful of GRBs in the prompt phase, and from GRB 
940217 in the afterglow phase. 

The photons received during the prompt phase are emit- 
ted closest to the central source; therefore a correct theory 
of the origin of the prompt emission reveals properties of the 
central engine. We contrast an internal shock/colliding shell 
model with an external shock model for the prompt 7-ray 
emission. Beaming breaks in afterglow light curves can be 
used to deduce the collimation angle, and hence the true en- 
ergetics of the GRB source. The most popular current model 
for the central sources of long-duration GRBs is the hyper- 
nova/collapsar model, where an evolved, massive star's cen- 
tral core collapses to a black hole while accreting material 
at the rate of w 0.1- a few Solar masses per second. The 
coalescence of neutron stars and black holes has been less 
successful as a model for the long-duration GRBs because 
host galaxy evidence indicate that a massive star origin is 
more likely; a compact-object coalescence model may still 
apply to the short duration GRBs. The supranova model of 
Vietri and Stella invokes a two-step collapse process where 
the core collapse of a massive star forms a rotationally stabi- 
lized neutron star following a supernova explosion. The neu- 
tron star later collapses to form a black hole. This model has 
several advantages over the hypernova/collapsar model to ac- 
count for the evidence of beaming and the X-ray absorption 
and line observations, including the presence of a highly en- 
riched circumburster medium (CBM), and a uniform energy 
scale set by the mass of the neutron star. 

Some important issues in GRB cosmology are mentioned, 
though the remainder of this review focuses on the role that 
GRBs may play in cosmic -ray origin studies. GRB sources 
release energy into an L* galaxy such as the Milky Way 
at the time averaged rate of ~ 10 40±1 ergs s^ 1 , which is 
nearly equal to the cosmic ray power into our Galaxy. Rel- 
ativistic blast waves can accelerate particles well beyond the 
knee of the cosmic-ray spectrum through second-order pro- 
cesses, overcoming difficulties that first-order Fermi acceler- 
ation faces to accelerate cosmic rays to energies above the 
knee of the cosmic ray spectrum in SNe. If the beaming re- 
sults are correct, then GRBs occur in our Galaxy at the rate of 
once every several thousand years, suggesting that GRBs are 
associated with a subset of SNe. Finally, we consider high- 
energy neutrino production from GRBs, and show how high- 
energy neutrino and gamma-ray observations of GRBs can 



discriminate between an internal and external shock model 
for the prompt gamma-ray emission, with important implica- 
tions for central engine theory. 

2 Observations and Phenomenology 

Because the emphasis here is on theory, we can only briefly 
summarize the basic features of GRB observations and phe- 
nomenology (for a recent review of the observations, see van 
Paradijs et al. (2000)). GRBs are brief flashes of radiation at 
hard X-ray and soft 7-ray energies that display a wide variety 
of time histories, though in <~ 25% of the cases a characteris- 
tic single-pulse profile is observed, consisting of a rapid rise 
followed by a slower decay (Fishman and Meegan, 1995). 
GRBs were first detected at soft 7-ray energies with wide 
field-of-view instruments. Peak soft 7-ray fluxes reach hun- 
dreds of photons cm~ 2 s _1 in rare cases. The BATSE in- 
strument is sensitive in the 50-300 keV band, and provides 
the most extensive data base of GRB observations during the 
prompt phase. It searches for GRBs by examining strings of 
data for > 5.5cr enhancements above background on the 64 
ms, 256 ms, and 1024 ms time scales, and triggers on GRBs 
as faint as w 0.5 ph cm~ 2 s _1 , corresponding to energy flux 
sensitivities £ 10~ 7 ergscm~ 2 s _1 . 

The integral size distribution of BATSE GRBs in terms of 
peak flux <fi p is very flat below <~ 3 ph cm~ 2 s _1 , and be- 
comes steeper than the —3/2 behavior expected from a Eu- 
clidean distribution of sources at <f> p >, 10 ph cm~ 2 s _1 (van 
Paradijs et al., 2000). GRBs typically show a very hard spec- 
trum in the hard X-ray to soft 7-ray regime, with a photon in- 
dex breaking from w — 1 at photon energies E p h ^ 50 ke V to 
a —2 to —3 spectrum at E p h £ several hundred keV (Band 
et al., 1993). Consequently, the distribution of the peak pho- 
ton energies E p k of the time-averaged vF v spectra of BATSE 
GRBs are typically found in the 100 keV - several MeV range 
(Mallozzi et al., 1995). The differential peak-flux size distri- 
bution (Meegan et al., 1996) and the E p k distribution (Mal- 
lozzi et al., 1997) are shown by the data points in the top and 
bottom panels of Fig. 1, respectively. 

The duration of a GRB is defined by the time during which 
the middle 50% (£50) or 90% (igo) of the counts above back- 
ground are measured. A bimodal duration distribution is 
measured, irrespective of whether the £50 or igo durations are 
considered (Kouveliotou et al., 1993). About two-thirds of 
BATSE GRBs are long-duration GRBs with t 90 > 2 s, with 
the remainder comprising the short-duration GRBs. The top 
panel of Fig. 2 shows the i 50 duration distribution of BATSE 
GRBs. 

The Beppo-SAX GRB observations revealed that essen- 
tially all long-duration GRBs have fading X-ray afterglows 
(Costa et al., 1999). The Wide Field Camera on Beppo-SAX 
has sensitivity down to <~ 10~ 10 ergs cm~ 2 s _1 with ^ 10' 
error boxes. Spacecraft slewing requires 6-8 hours, but per- 
mits Narrow Field Instrument X-ray observations with sen- 
sitivity down to ~ 10~ 14 ergs cm~ 2 s _1 and error boxes 
^,0.5'. The first X-ray afterglow was obtained from GRB 
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Fig. 1. Differential distribution of the number of GRBs with peak 
fluxes measured on the 1024 ms time scale with BATSE is plotted 
in the upper panel, and the distribution of the photon energies of the 
vF v peaks are plotted on the bottom panel. Histogram gives a fit to 
the data using an external-shock model for the prompt phase. 



Fig. 2. Data points in the top panel give the t$o duration distribution. 
External shock-model fit to the long duration distribution is given 
by the histogram. Bottom panel shows the redshift distribution of 
BATSE GRBs predicted by the model (Bottcher and Dermer, 2000). 



970228 (Costa et al., 1997), which revealed an X-ray source 
which decayed according to a power-law behavior <f>x oc t x , 
with x ~ —1.33. Typically, \ ~ —1.1 to —1.5 in X-ray 
afterglow spectra. 

The small X-ray error boxes allow deep optical and ra- 
dio follow-up studies. GRB 970228 was the first GRB from 
which an optical counterpart was observed (van Paradijs et 
al., 1997), and GRB 970508 was the first GRB for which 
a redshift was measured (Djorgovski et al., 1997; Metzger 
et al., 1997). Detection of optical emission lines from the 
host galaxy, and absorption lines in the fading optical after- 
glow due to the presence of intervening gas has provided red- 
shifts for about 20 GRBs (Djorgovski et al., 2001). No opti- 
cal counterparts are detected from approximately one-half of 
GRBs with well-localized X-ray afterglows, and are termed 
"dark" bursts. These sources may be undetected in the op- 
tical band because of dusty media (Groot et al., 1998). Ap- 
proximately 40% of GRBs have radio counterparts, and the 
transition from a scintillating to smooth behavior in the radio 
afterglow of GRB 980425 provides evidence for an expand- 
ing source (Frail et al., 1997). Fig. 3 gives the distribution of 
redshift and apparent isotropic 7-ray energy releases for 17 
GRB for which the redshift is known, using the compilation 
of Frail etal. (2001). 



X-ray emission lines, possibly due to Fe Ka fluorescence, 
were detected during a rebrightening phase in the afterglow 
of GRB 970508 (Piro et al., 1999), and in the afterglow spec- 
tra of GRB 991216 (Piro et al., 2000) and GRB 000214 (An- 
tonelli et al., 2000). A marginal detection of an X-ray emis- 
sion feature was also reported from GRB 970828 (Yoshida 
et al., 1999). A transient absorption feature was detected 
from GRB 990705 (Amati et al., 2000), and X-ray absorp- 
tion in excess of the Galactic hydrogen column density has 
also been reported in GRB 980329 (Frontera et al., 2000). 
These results provide an alternative method of redshift deter- 
mination, and provide important clues about the environment 
of the progenitor object, showing that large quantities of iron 
must be present in the vicinity of these sources. 

The EGRET instrument on the Compton Observatory de- 
tected 7 GRBs with > 30 MeV emission during the prompt 
phase (Catelli et al., 1998), including the extraordinary burst 
GRB 940217, which displayed > 100 MeV emission 90 min- 
utes after the onset of the GRB, including one photon with 
energy near 20 GeV (Hurley et al., 1994). This gives unam- 
biguous evidence for the importance of nonthermal processes 
in GRBs. TeV radiation has been reported to be detected with 
the Milagrito water Cherenkov detector from GRB 970417a 
(Atkins et al., 2000). If correct, this requires that this source 
be located at z ^ 0.3 in order that 77 attenuation with the 
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Fig. 3. Distribution of redshifts and apparent isotropic 7-ray energy 
releases. 



diffuse intergalactic infrared radiation to be small. 

A class of X-ray rich GRBs, with durations of order of 
minutes and X-ray fluxes in the range 10~ 8 -10~ 7 ergs cm -2 
s _1 in the 2-25 keV band, has been detected with many X- 
ray satellites, including Ariel V, HEAO-1, ROSAT, Ginga, 
and Beppo-SAX (Heise et al., 2001). We also note sev- 
eral phenomenological correlations, including claims of an 
inverse correlation of luminosity and BATSE channel lags 
(Norris et al., 2000), correlations between the duration of 
quiescent and subsequent pulse periods in separated pulses 
(Ramirez-Ruiz et al., 2001), and indications that the peak lu- 
minosity is correlated with variability (Reichart et al., 2001). 
Other phenomenological behaviors in the prompt emission 
include hard-to-soft evolutions and hardness-intensity corre- 
lations, and a quantitative relation between integrated fluence 
and E p k in well-defined pulses (Liang and Kargatis, 1996). 



3 Elementary Blast Wave Theory 

The fireball/blast-wave model for GRBs was proposed by 
Rees and Meszaros (1992) and Paczyhski and Rhoads (1993), 
and successfully predicted power-law afterglows (Meszaros 
and Rees, 1997). It has emerged as the standard model to 
interpret GRB prompt and afterglow emissions (for a recent 
review, see Meszaros (2002)). The simplest version of the 
standard blast wave model involves a spherical, uncollimated 
explosion taking place in a uniform surrounding medium. A 
relativistic pair fireball is formed when an explosion deposits 
a large amount of energy into a compact volume. The pres- 
sure of the explosion causes the fireball to expand, with the 
thermal kinetic energy of the explosion being transformed 
into directed kinetic energy due to strong adiabatic cooling 
of particles in the comoving frame (Meszaros et al., 1993; 
Piran, 1999). Because of the Thomson coupling between 
the particles and photons, most of the original explosion en- 



ergy is carried by the baryons that were originally mixed 
into the explosion. Under certain conditions involving less 
energetic, temporally extended, or very baryon-clean explo- 
sions, neutrons can decouple from the flow (Derishev et al., 
1999). If this does not occur, then the coasting Lorentz fac- 
tor To = Eo/Mi>c 2 , where M& is the baryonic mass and 
Eo = 10 52 i?52 ergs is the apparent isotropic energy release. 

In the simplest version of the model, the blast wave is ap- 
proximated by a uniform thin shell, and particle accelera- 
tion is assumed to take place at the forward shock only. A 
forward shock is formed when the expanding shell accel- 
erates the external medium, and a reverse shock is formed 
due to deceleration of the cold shell. The forward and re- 
verse shocked fluids are separated by a contact discontinuity 
and have equal kinetic energy densities. From the relativistic 
shockjump conditions (Sari and Piran, 1995), 4r(T— l)no = 
4f (f — l)n' sh , where T is the blast wave Lorentz factor, f 
is the Lorentz factor of the reverse shock in the rest frame 
of the shell, n' sh is the density of the unshocked fluid in the 
proper frame of the expanding shell, and tiq is the density of 
the CBM, here assumed to be composed of hydrogen. Parti- 
cle acceleration at the reverse shock is unimportant when the 
reverse shock is nonrelativistic, which occurs when 




(1) 



The unique feature of GRBs is that the coasting Lorentz 
factor Tq may reach values from hundreds to thousands. By 
contrast, Type la and II SNe have ejecta speeds that are in the 
range of ~ 5000-30000 km s _1 . The relativistic motion of 
the radiating particles introduce many interesting effects in 
GRB emissions that must be properly taken into account. 

Three frames of reference are considered when discussing 
the emission from systems moving with relativistic speeds: 
the stationary frame, which are denoted here by asterisks, 
the comoving frame, denoted by primes, and the observer 
frame. The differential distance traveled by the expanding 
blast wave during differential time dt* is simply dx — f3cdt*, 
where (3 = yl — Y~ 2 . Due to time dilation, dx = fSYcdt' . 
Because of time dilation, the Doppler effect, and the cosmo- 
logical redshift z, the relationship between comoving and ob- 
server times is (l + z)Tdt'(l- (3cos6) = (l + z)dt'/5 = dt, 
where 9 is the angle between the emitting element and the ob- 
server and 5 — [T(l — (3 cos 6>)] _1 is the Doppler factor. For 
on-axis emission from a highly relativistic emitting region, 
we therefore see that dt = (1 + z)dx/T 2 c; consequently the 
blast wave can travel a large distance T 2 cAt during a small 
observing time interval. A photon measured with dimension- 
less energy e = hv/m e c 2 is emitted with energy 5e' j (1 + z). 

A blast wave expanding into the surrounding medium acts 
as a fluid if there is a magnetic field in the comoving frame 
sufficiently strong to confine the particles within the width of 
the blast wave. The blast-wave width depends on the dura- 
tion of the explosion and the spreading of the blast-wave par- 
ticles. Because sub-second structure is regularly observed 
in GRB light curves, the central engine must be ;$ 1 lt-s, 
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and is probably much thinner in view of 1-10 ms variabil- 
ity observed in some GRBs (Bhat et al., 1992). We denote 
this length scale by Ao cm. The shell will spread radially 
in the comoving frame by an amount x'n ~ v spr t', where 
t' = x/((3Tc) is the available comoving time and v spr is 
the spreading speed. This implies a shell width in the ob- 
server frame of A ~ x/Tq due to length contraction, and 
a spreading radius x spr = TqA , assuming that the shell 
spreads with speed v spr = c. The width of the unshocked 
blast-wave fluid in the rest frame of the explosion is there- 
fore A = min(A , x/Tg) (Meszaros et al., 1993; Sari et al., 
1996). Milligauss fields are sufficient to confine relativistic 
electrons. 

As the blast wave expands, it sweeps up material from the 
surrounding medium to form an external shock (Meszaros 
and Rees, 1993). In a colliding wind scenario, the blast wave 
intercepts other portions of the relativistic wind (Rees and 
Meszaros, 1994). Protons captured by the expanding blast 
wave from the CBM will have total energy Tm p c 2 in the 
fluid frame, where m p is the proton mass. The kinetic en- 
ergy swept into the comoving frame by an uncollimated blast 
wave at the forward shock per unit time is given by 



dE' 

~d¥ 



= 4irx 2 n a m p c 3 f3r(T - 1) 



(2) 



(Blandford and McKee, 1976). One factor of T represents the 
increase of external medium density due to length contrac- 
tion, the factor (T — 1) is proportional to the kinetic energy of 
the swept up particles, and the factor (3 is proportional to the 
rate at which the particle energy is swept. Thus the power is 
oc T 2 for relativistic blast waves, and oc /3 3 for nonrelativistic 
blast waves. This process provides internal energy available 
to be dissipated in the blast wave. 

A proton that is captured by and isotropized in the comov- 
ing frame will have energy Tm p c 2 in the comoving fluid 
frame, or energy T 2 m p c 2 in the observer frame. The ex- 
panding shell will therefore begin to decelerate when Eq = 
T M b c 2 = rgm p c 2 (47ra^no/3), giving the deceleration ra- 
dius 

X d EE )V3 a 2 . 6 x l 16 ( ^_)l/3 cm (3) 

(Rees and Meszaros, 1992; Meszaros and Rees, 1993), where 
T300 = ro/300. Acceleration at the shock front can in- 
ject power-law distributions of particles. In the process of 
isotropizing the captured particles, magnetic turbulence is 
introduced (Pohl and Schlickeiser, 2000) that can also accel- 
erate particles to very high energies through a second-order 
Fermi process. 

The deceleration time as measured by an on-axis observer 
is given by 



t d =(l + z) 



Xd 



9.6 (1 + z) E x 



52 



/? I> 



A) 



Tf, 



1/3 



300 n 



(4) 



The factor [3 Q = 1/ y 1 — T generalizes (Dermer and 
Humi, 200 1 ) the result of Meszaros and Rees ( 1 993) for mildly 



relativistic and nonrelativistic ejecta, as in the case of Type 
la and Type II supernova explosions. The Sedov radius 



is = T 2 fx d = ( 



3£ 



47rra p c 2 no ' 



vl/3 ca 
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1.2 x 10 18 (^)V3 cm 6.6 x lO 18 ^) 1 / 3 cm , 
n n 

where m Q is the total (rest mass plus kinetic) explosion en- 
ergy in units of Solar rest mass energy. For relativistic ejecta, 
is refers to the radius where the blast wave slows to mildly 
relativistic speeds, i.e., r <~ 2. The Sedov radius of a SN that 
ejects a 10 Mq envelope could reach <~ 5 pc or more. The 
Sedov age t s = £ s /v = 700(m Q /n ) 1/3 /(uo/0.01c) yrfor 
nonrelativistic ejecta, and is equivalent to td in general. 

The evolution of an adiabatic blast wave in a uniform sur- 
rounding medium for an explosion with a nonrelativistic re- 
verse shock is given by (Dermer and Humi, 2001) 



P(x) = 



Po 



A>ro j 



x <C x d 



y/1 + {X/Xdf l /3o (x ) -3/2 5 Xd< , x 



, (6) 



where P — (3T and P — Po^o represent dimensionless bulk 
momenta of the shocked fluid. Equation (6) reduces to the 
adiabatic Sedov behavior for nonrelativistic (Po^o -C 1) ex- 
plosions, giving v oc x~ 3 l 2 , x oc i 2 / 5 , and v oc i -3 / 5 , as is 
well-known. In the relativistic (To ^> 1) limit, T oc x~ z / 2 
and t = c _1 J dx/T 2 oc J dx x 3 , yielding x oc i 1 / 4 and T oc 
t^ 3 / 8 when r» 1. The adiabatic behavior does not apply 
when radiative losses are important. A nonrelativistic super- 
nova shock becomes radiative at late stages of the blast-wave 
evolution. It is not clear whether a GRB fireball is highly ra- 
diative or nearly adiabatic during either its gamma-ray lumi- 
nous prompt phase or during the afterglow (Waxman, 1997; 
Vietri, 1997). In the fully radiative limit, equation of energy 
conservation gives Po^oMq — /?r[M + 47rno:r 3 /3], giving 
the asymptotes T oc x~ 3 when r > T » 1 and (3 oc x~ 3 
when r - 1< 1. 

Most treatments employing blast-wave theory to explain 
the observed radiations from GRBs assume that the radiating 
particles are electrons. The problem here is that <~ m p /m e <~ 
2000 of the nonthermal particle energy swept into the blast- 
wave shock is in the form of protons or ions unless the sur- 
roundings are composed primarily of electron-positron pairs 
(Konigl and Granot, 2001). For a radiatively efficient system, 
physical processes must therefore transfer a large fraction of 
the swept-up energy to the electron component. In elemen- 
tary treatments of the blast-wave model, it is simply assumed 
that a fraction e e of the forward-shock power is transferred 
to the electrons, so that 



dE' 



(7) 



If all the swept-up electrons are accelerated, then joint nor- 
malization to power and number gives 



7™n = e e (^4)(^)(r-l) =e e k p (^)T 
p — 1 m e m e 



(8) 
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for 2 < p < 3, where the last expression holds when T >• 1 
andfc p = (f>-2)/(p-l). 

The strength of the magnetic field is another major uncer- 
tainty. The standard prescription is to assume that the mag- 
netic field energy density Uf, = B 2 /8tt is a fixed fraction 
e B of the downstream energy density of the shocked fluid. 
Hence 



localized emission region, and the received vF v spectrum is 
given by 
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B 2 

— = 4e B n m p c 2 (T 2 - T) . 

87T 



(9) 



It is also generally supposed in simple blast-wave model cal- 
culations that some mechanism - probably the first-order shock 
Fermi process - injects electrons with a power-law distribu- 
tion between electron Lorentz factors 7 m i„ < 7 < 72 down- 
stream of the shock front. The electron injection spectrum in 
the comoving frame is modeled by the expression 



dt , d ^ =K el p , for lmm < 7 < 72 , 



(10) 



where K e is normalized to the rate at which electrons are 
captured. 

The maximum injection energy is obtained by balancing 
synchrotron losses and an acceleration rate given in terms of 
the inverse of the Larmor time scale through a parameter e2, 
giving 



Here the observed flux is proportional to four powers of the 
Doppler factor: two associated with solid angle, one with 
energy and one with time. 

From this formalism, analytic and numerical models of 
relativistic blast waves can be constructed. It is useful to 
define two regimes depending on whether j m i n < 7c, which 
is called the weak cooling regime, or "f m i n > -f c , called the 
strong cooling regime (Sari et al., 1998). If the parameters p, 
e e and e B remain constant during the evolution of the blast 
wave, then a system originally in the weak cooling regime 
will always remain in the weak cooling regime. In contrast, 
a system in the strong cooling regime will evolve to the weak 
cooling regime. 

For a power-law injection spectrum given by equation (10), 
the cooling comoving nonthermal electron spectrum can be 
approximated by 



Kb) 



N° 70 s ' 1 h * ' 7o ^ 7 ^ 7i 

l7f +1 - S 7- (p+1) , 71^7^72, 



(15) 



72 S 4 x 10 7 e 2 /v/B(G) . 



(11) 



A break is formed in the electron spectrum at cooling elec- 
tron Lorentz factor -f c , which is found by balancing the syn- 
chrotron loss time scale t' syn with the adiabatic expansion 

giving 



and N® = Awx 3 no/3 for the assumed system. In the weak 
cooling regime, s = p, 70 = jmin and 71 = j c , whereas in 
the strong cooling regime, s = 2, 70 = j c and 71 = ■y m in- 

As an example of the elementary theory, consider the tem- 
poral index in the strongly cooled regime for high energy 
electrons with 7 > 71 = -f m oc V, 70 = j c oc l/(T 3 t), and 



7c = 



3m e 



16e Bnom p caT^ 3 t 



(12) 



(Sari et al., 1998). For an adiabatic blast wave, T oc t~ 3 / 8 , 
so that j m i n oc i~ 3 / 8 and 7 C oc i 1 / 8 . 

The observed vF v synchrotron spectrum from a GRB de- 
pends on the geometry of the outflow. Denoting the co- 
moving spectral luminosity by L' syn (e') = e'(dN' /de'dt'), 
then e'L' syn (e') = ^u B c(7T7 3 N e (-/), with 7 = yje' /e B and 
e B = B/B cr = 5/(4.41 x 10 13 G). For a spherical blast- 
wave geometry, the spectral power is amplified by two pow- 
ers of the Doppler factor <5 for the transformed energy and 
time. The vF v synchrotron spectrum is therefore 



2F 2 



{4c<7 T B 2 ~ fc /24irm e c 2 )- 1 , s = 2. From equation (13), f e oc T 4 x 3 7^ 1 7[ +1 ~M 2 -p)/ 2 

/r 2 ~P oc T 2{ P- lS >x 3 e {2 -^/ 2 /t oc t x e*, with x = (2 - 3p) /4 
and a„ = (2 — p)/2. A few other such results include 
the decay of a strong and weak cooling vF v peak frequency 
e p k oc i~ 3 / 2 and e p k oc i~ 3p / 2 , resepctively. A cooling in- 
dex change by A(a„) = 1/2 is accompanied by a change 
of temporal index from \ = 3(1 — p)/4 to x = (2 — 3p)/4 
at late times, so that Ax = 1/4. Extension to power-law 
radial electron profiles is obvious, and beaming breaks intro- 
duce two additional factors of T from S m T in equation (14) 
when the observer sees beyond the causally connected re- 
gions of the jetted blob, noting that N® must be renomalized 
appropriately. 

Fig. 4 shows spectral and temporal indices that are derived 
from the preceding analytic considerations of a spherical adi- 
abatic blast wave that is in the strong cooling regime during 
a portion of its evolution (Dermer and Bottcher, 2002). We 
use the notation that the vF v spectrum is described by 



(l + z)e 
2Te B 



(13) 



fe ex t*e c 



(16) 



where c£l is the luminosity distance. 

For collimated blast waves with jet opening angle 6j, equa- 
tion (13) does not apply when T < 1/6 j because portions of 
the blast wave's radiating surface no longer contribute to the 
observed emission. In this case, the blast-wave geometry is a 



Observations at a specific photon energy e will detect the 
system evolving through regimes with different spectral and 
temporal behaviors. The relativistic Sedov phase corresponds 
to the afterglow regime, and we also consider in this figure a 
possible relativistic reverse shock (RRS) phase. If equation 



7 



log e 



I \ 


k-3 

<x r vi 

X = (2-p)/2 


Strong Cooling 

.... 


\ oct° 
3 \ 2 


(3-p)/2 \l-(p/2) 

t° \ 
oc t \ 

11/3 \2 


a = (2-p)/2 

I- * 

1/2 ~" - min 


(2-3p)/4 


4/3 \l /2 
3 \ 


1/2 ^ 


(2-p)/2 

\ 

-1/4 \ 


4/3 


e 
"""•»£. 

4/3 - 


1/2 \ 

\ t -3/2 3(1-p)/4 
f 1 (3-p)/2 
\ / 




| 4/3 ^ 


J v / > 

| 1/2 \^~~~-^> 


Prompt phase 


1/6 
4/3 

RRS phase 
I I 


\ K t 
\ 

Relativistic Sedov phase 
> 



logt NR log t d logt 4 log t 



Fig. 4. Possible time profiles and spectral behavior of the nonther- 
mal synchrotron radiation that could arise for a system that is in the 
fast cooling regime during a portion of the adiabatic evolution of a 
blast wave in a uniform surrounding medium. The temporal indices 
X and vF v spectral indices ct v are shown, respectively, above and 
below the lines that represent different families of possible emission 
trajectories. 



(1) is satisfied with T replaced by To, then the blast wave 
will not evolve through the RRS phase. Even for this simple 
system, a wide variety of behaviors are possible. Inclusion of 
additional effects, including a nonuniform external medium 
(Meszaios et al., 1998) or blast-wave evolution in a partially 
radiative phase (Bottcher and Dermer, 2000a) will introduce 
additional complications. 

The heavy solid and dotted lines in Fig. 4 represent the 
evolution of the cooling photon energies e c and e m i n , respec- 
tively. In the lower-right hand corner of the diagram, one sees 
that observations in the afterglow phase may detect a transi- 
tion from the uncooled portion of the synchrotron spectrum 
where x = (3 — p)/2 and a u = (3 — p)/2 to a cooled por- 
tion of the synchrotron spectrum where \ — (2 — 3p)/4 and 
a u = (2 — p)/2. Thus a change in photon index by 0.5 units 
is accompanied by a change in temporal index by 1/4 unit, in- 
dependent of p. Comparison with afterglow observations can 
be used to test the blast-wave model (Galama et al., 1998). 
The upcoming Swift telescope will be able to monitor be- 
haviors in the time interval between the prompt and afterglow 
phase, and to search for other regimes of blast wave evolution 
and evidence for evolution in the non-adiabatic regime. 

Fig. 5 shows a numerical model of a relativistic blast wave 
calculated with the code developed by Chiang and Dermer 
(1999), which takes into account light travel-time delays from 
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Fig. 5. Calculations of spectral energy distributions emitted by a 
relativistic blast wave with initial Lorentz factor To = 600 that 
is energized by sweeping up material from a uniform surrounding 
medium with density n. Curves are labeled by the base 10 logarithm 
of the observing time in seconds. The analytic and numerical mod- 
els are shown, with the latter curves identified by data points. Other 
parameters of the calculation are shown in the legend, based on pa- 
rameters taken from the fit to GRB 970508 by Wijers and Galama 
(1999). Upper panel shows spectral energy distribution in a vL v 
representation, and the lower panel shows light curves at different 
energies in an L v representation. 



different portions of the blast wave, synchrotron and syn- 
chrotron self-Compton emission, as well as self-absorption. 
The numerical results are compared with analytic spectral 
energy distributions (Dermer et al., 2000) in which an SSC 
component is included in the approach of Sari et al. (1998). 
We use parameters derived by Wijers and Galama (1999) to 
fit afterglow data from GRB 970508. Here the total energy 
injected into the fireball is E — 3.5 x 10 52 ergs, and the 
density of the external medium is n — 0.03 cm -3 . The SSC 
component becomes increasingly important with increasing 
e e /e B (Sari et al., 1996; Moderski et al., 2000), and will 
be observed as features in the X-ray and optical light curves 
(Dermer et al., 2000a; Zhang and Meszaios, 2001). As can 
be seen, the analytic expressions provide a reasonable repre- 
sentation of the numerical model in the optical/X-ray regime, 
but the discrepancies at radio frequencies are significant, as is 
clear from the second panel in Fig. 5. Although the times of 
the temporal breaks arising from synchrotron emission emit- 
ted by the lowest energy electrons are in good agreement 
with the numerical calculation, the analytic model can be dis- 



8 



2 10 



1.6 10 



V 1.2 10° - 



8 10 



4 10 



10" 




10 20 30 40 50 60 70 80 
t (sec) 

Fig. 6. Numerical calculations of light curves measured at different 
observed photon energies. These profiles correspond to the most 
common prompt GRB light curves. The hard-to-soft evolution is 
indicated by the variation in the photon index. For this model, an 
explosion with total energy Eq = 10 54 ergs and coasting Lorentz 
factor To = 300 expands into a uniform surrounding medium with 
density no = 100 cm 2 . 




(a) 



r = 1000 

o 

(relative flux divided by 5) 



r o =ioo 



40 60 

time (s) 



80 



100 



Fig. 7. Light curves calculated at 100 keV from a fireball with other 
parameters the same as in Fig. 6. 



crepant in flux values by an order-of-magnitude or more, so 
that any quantitative statement based upon analytic models 
must be made with care, particularly in the radio regime. 

Fig. 6 shows light curves calculated from an adiabatic blast 
wave that expands into a uniform surrounding medium (Der- 
mer et al., 1999b). The light curves display a temporal pro- 
file that rapidly rises and slowly decays, and reproduces the 
so-called fast-rise, exponential decay profiles observed in a 
significant fraction of GRBs. This indicates that the prompt 
emission in these GRBs, at least, arises from an external 
shock model. An external shock origin also accounts for the 
basic phenomenology of the GRBs with smooth profiles, in- 
cluding the hard-to-soft evolution and the hardness-intensity 
correlation observed in GRB light curves. 

3.1 Clean and Dirty Fireballs 

If a subset of prompt GRB light curves originate from an ex- 
ternal shock model in the prompt phase, one might ask what 
this means for the remainder of the GRBs. A majority of 
the GRB community believes that an active central engine 
is required to reproduce the spiky, short timescale variablity 
observed in a large number of GRB light curves, and that an 
impulsive explosion cannot produce the temporal gaps de- 
tected in some GRB light curves. Inhomogeneities in the 
surrounding medium can, in fact, produce short timescale 
variability and gaps in the light curves, as demonstrated by 
Dermer and Mitman (1999). More serious issues for the ex- 
ternal shock model for the prompt emission are posed by the 
observations of optical emission observed in GRB 990123 in 
the prompt phase if due to reverse shock emission (Akerlof et 
al., 1999), and the correlated durations of quiescent and sub- 
sequent peak durations in GRB light curves (Ramirez-Ruiz 
et al., 2001). 

Rather than revisit this controversy, it seems more useful 



to focus on the predictive capability of the external shock 
model. Because of the small number of parameters, it is sim- 
ple to outline the generic behavior of GRBs in terms of the 
model's most sensitive parameters, namely the total explo- 
sion energy Eq, the surrounding medium density no, and the 
coasting Lorentz factor r . 

The basic features of GRB fireballs can be summarized by 
a few simple equations (Dermer et al., 1999a). The first is 
that the total energy liberated in photons during the prompt 
phase is 

E 1 = 0.1 U t d , (17) 

where the precise value of the coefficient, here taken to be 
10%, depends on how radiative the fireball is during this 
phase. The power n during the prompt phase is 

n cx (n £ 2 ) 1/3 if 3 (18) 
and, as already seen in equation (4), the deceleration time 

t d oc(^)V3r 8 / 3 . (19) 

Most of the emission during the prompt phase is radiated at 
photon energy 

E pk oc nJ /2 ro . (20) 

The strong dependence of these quantities on To is evi- 
dent. A detector observing in a fixed range of photon en- 
ergy centered at Edet, such as BATSE, will have strong se- 
lection biases to detect fireballs with different r factors. The 
biases against detecting dirty fireballs, namely those with 
small values of r and large baryon-loading, arise because 
the peak luminosity measured by a detector declines rapidly 
when r f . The quantity f is the Lorentz factor of a 
fireball that produces the bulk of its emission near the sen- 
sitive range of a detector, and is defined by E p k = Edet oc 
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rip 2 f g, where E det is the photon energy to which the GRB 
detector is largest effective area. The peak luminosity is 
larger for clean fireballs with T > f and small baryon 
loading, but the shorter durations of the clean events mean 
that detectors trigger on fiuence rather than peak flux when 
r ^> f . GRB detectors are therefore also biased against 
the detection of clean fireballs. This explains the tendency of 
detectors to detect GRBs that produce prompt radiation with 
vF v peak in the waveband at which the detector has greatest 
sensitivity. 

These considerations resolve the apparent paradox that a 
beamed scenario should seem to produce a wide range of 
E pk , contrary to observations (Dermer et al., 1999b) (see 
lower panel in Fig. 1). To avoid fine-tuning of Tq, one re- 
alizes that there must be two additional classes of explosive 
phenomena that remain to be discovered by detectors with 
appropriate design (Dermer et al., 1999a), namely the large 
baryon-load, low Tq "dirty" fireballs, and the low baryon- 
load, high T "clean" fireballs. 

This is apparent from Fig. 7, which shows how a variation 
of a factor of 3 in T from 300 to 100 can make a GRB es- 
sentially undetectable in the BATSE band. The low Lorentz 
factor, dirty fireballs are more extended in time, and radiate 
their energy at lower observed photon frequencies. These 
may correspond to the X-ray rich GRBs that have been de- 
tected with Beppo-SAX and previous X-ray detectors (Heise 
et al., 2001). 

By taking the triggering properties of BATSE into account, 
the external shock model can be used to fit the statistical 
data on GRBs (Bottcher and Dermer, 2000). The results are 
shown in Figs. 1 and 2, where simple power-law distributions 
of Tq and Eq were assumed. The underlying assumption of 
this model is that the rate of GRBs per comoving density fol- 
lows the star formation history of the universe. These results 
show that fine-tuning is not required to explain the tendency 
of Epk to be within a narrow range. The tendency of softer 
bursts to have longer durations is a direct consequence of this 
model, though only the long duration GRBs are fit in this 
model. A wide range of apparent isotropic energy releases 
Eq is used to fit the statistical data. If GRBs are beamed, the 
distribution of Eq would be consistent with recent analyses 
that indicates that GRBs have a standard energy reservoir, 
with the apparent variation of E reflecting the range of jet 
opening angles. 

3.2 Beaming in GRBs 

An observer will receive most emission from those portions 
of a GRB blast wave that are within an angle ~ l/T to the 
direction to the observer. As the blast wave decelerates by 
sweeping up material from the CBM, a break in the light 
curve will occur when the jet opening half-angle 9j becomes 
smaller than l/T. This is due to a change from a spherical 
blast wave geometry, given by equation (13), to a geometry 
defined by a localized emission region, as given by equation 
(14). Assuming that the blast wave decelerates adiabatically 
in a uniform surrounding medium, the condition 9j = l/T = 
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Fig. 8. Light curves calculated at various observing energies and in- 
clination angles 6 for a fireball with parameters the same as in Fig. 6 
and jet opening half-angle 8j = 10°. The dashed, dotted, and solid 
curves give light curves measured at 3 keV, V band, and 8.6 GHz 
radio frequencies, and the observing angles 9 = 0°, 10°, 12°, 30°, 
and 90°. 
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can be derived (Sari et al., 1999). Note that the beaming 
angle is only weakly dependent on no and Eq. 

Fig. 8 shows calculations of a GRB blast wave in a beamed 
geometry (Dermer et al., 2000a) with 9j = 10° at different 
observing angles and energies. The jet emission is assumed 
to be uniform across its surface and not to spread laterally 
(see Rhoads (1999) when this is not the case). For the pa- 
rameters used here (£52 = 100, e# = 10~ 4 , e e = 0.5, 
«o = 100 cm~ 3 , p = 2.5), the beaming break occurs at 
~ 10 5 s, but is obscured at X-ray and optical frequencies by 
the appearance of an SSC component. The appearance of 
clearly defined breaks in the light curves of some GRBs, if 
due to beaming rather than to a variation in the density of the 
surrounding medium, limits the ratio e e /es that defines the 
strength of the SSC component (Dermer et al., 2000a). So- 
called "orphan" afterglows, which are those GRBs with jets 
pointed away from our line-of-sight which become visible 
upon deceleration of the blast wave, may be confused with 
dirty fireballs (Dermer et al., 2000a; Huang, et al., 2002). 

Afterglow modeling of multiwavelength spectra of GRB 
980703, GRB 990123, GRB 9905010, and GRB 991216 by 
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Panaitescu and Kumar (2001) shows that the CBM is more 
consistent with uniform surroundings than with a wind pro- 
file with density n(r) oc r~ 2 . They derive magnetic field 
parameters in the range 10~ 4 £ es ^ 0.05, electron energy- 
transfer parameters in the range 0.01 ;$ e e <, 0.1, and obtain 
low density surroundings (<~ 10~ 4 ;$ n <, 10 cm -3 ) in their 
fits. The jet angles Oj range from 1° to 4°. The beaming fac- 
tor for a one-sided jet is ss 1.3 x 10 4 / (9 ) 2 , so that many mis- 
aligned GRBs sources should exist for every detected GRB. 

Frail et al. (2001) have considered a larger sample of GRBs 
and show that the beaming results provide evidence for a 
standard energy reservoir. Letting r] 7 represent the efficiency 
for transforming total energy E to t (which should be distin- 
guished from the apparent isotropic energy release Eq) to 7 
ray energy, then E to t — ^-E^Ciso)/^^), where i? 7 (iso) is 
the apparent isotropic energy release in the form of GRBs. 
They find that the mean energy release by the GRB sources 
is (E GRB ) = 3x 10 51 /(?7 7 /0.2) ergs, and that the actual 
event rate of GRBs is s=a 500 x the observed rate. If correct, 
this result has important implications for the nature of the 
sources of GRBs. 



4 Source Models 

Keeping in mind that only members of the class of long- 
duration GRBs have measured redshifts, considerable evi- 
dence connects GRBs to star-forming regions (Djorgovski et 
al., 2001; Sokolov et al., 2001) and, consequently, to a mas- 
sive star origin. For example, the associated host galaxies 
have blue colors, consistent with galaxy types that are un- 
dergoing active star formation. GRB counterparts are found 
within the optical radii and central regions of the host galax- 
ies (Bloom, et al., 1999a). Lack of optical counterparts in 
some GRBs could be due to extreme reddening from large 
quantities of gas and dust in the host galaxy (Galama and Wi- 
jers, 2001). Supernova-like emissions have been detected in 
the late-time optical decay curves of a few GRBs (Reichart, 
1999; Bloom, et al., 1999b), and X-ray evidence for Fe K«- 
line signatures indicates that large quantities of highly en- 
riched material are near the sources of GRBs. The spatial 
and temporal coincidence of GRB 980425 with SN 1998bw, 
a Type Ic supernova, if true, also connects the sources of 
GRBs to SNe (Kulkarni et al., 1998; Pian et al., 1999). 

4. 1 Coalescing Compact Objects 

A central attraction of this model is that binary neutron stars, 
such as the Hulse-Taylor pulsar PSR 1913+16, exist. This 
system has a merger time scale w 3 x 10 8 yrs, which is much 
less than the Hubble time, so that these events undoubtedly 
occur. The timescale for the final coalescence event is very 
short, with the bulk of the energy released over a few tens 
of milliseconds (Janka et al., 1999). Such events should, 
however, occur in both spiral and elliptical galaxies, whereas 
the host galaxies of known GRBs generally exhibit active 
star formation and are therefore associated with spiral or dis- 



turbed galaxy hosts. The absence of counterparts far from the 
disks of galaxy hosts also conflicts with this model, as such 
old stellar systems can travel great distances before coales- 
cence. In particular, they should leave dusty regions where 
young and massive stars are born. 

Merger-rate estimates of neutron-star/neutron-star in addi- 
tion to neutron-star/black-hole binaries range from s=a 10~ 5 - 
10~ 6 yr~ 4 in our Galaxy, but are quite uncertain due to the 
sensitive dependence of merging timescale on initial binary 
separation (Narayan et al., 1991; Tutukov and Yungelson, 
1993). Even so, if GRBs outflows are highly beamed, event 
rates ^ 10~ 4 yr~ 4 are implied (see §6.3). Taken together, 
these various lines of evidence do not support a model of co- 
alescing compact objects for the long duration GRBs. 

4.2 Hypernova/Collapsar Model 

In the original formulation of this model (Woosley, 1993), a 
collapsar was a "failed" supernova in the sense that a core- 
collapse event failed to collapse to form a neutron star and 
instead produced a black hole. In order to produce a long- 
duration GRB with complex pulse structure, Woosley (1993) 
argued that accretion had to proceed over a period of time 
comparable to the prompt phase of a GRB. The progeni- 
tor star was suggested to be a rotating Wolf-Rayet star that 
produces, upon collapse, an accretion disk of several tenths 
of Solar masses. Because of the accretion geometry, jetted 
emission from neutrino annihilation was formed along the 
axis of the system. 

Paczyhski (1998) introduced the word "hypernova" to re- 
fer to the extremely luminous prompt GRB and afterglow 
events, and concurred with Woosley that GRBs are formed 
by the collapse of a massive star to a black hole surrounded 
by a massive disk and torus. He also suggested that the pre- 
cursor X-ray emission observed with Ginga in several GRBs 
(Murakami et al., 1991) was a signature of a dirty fireball 
produced as the ejecta cleared away overlying baryonic mate- 
rial, after which a highly relativistic, clean fireball could then 
emerge. In the version of the hypernova model advanced by 
Paczyhski (1998), the events are in fact as energetic as 10 54 
ergs and 10 4 -10 5 times rarer than SNe. 

Hydrodynamical simulations of collapsars have specifically 
treated, for example, the evolution of a 35 M Q main-sequence 
star whose 14 Mq helium core collapses to form a 2-3 Mq 
black hole (MacFadyen and Woosley, 1999). Provided that 
the core has a large amount of angular momentum, a de- 
layed accretion event can be formed by the infalling mat- 
ter. The black-hole jet simulation is performed by inject- 
ing <~ 10 51 ergs s _1 , which is argued to be generated by 
magneto-hydrodynamical processes rather than neutrino an- 
nihilation, insofar as the latter process produces inadequate 
energy to break through the core. In cases where break- 
through is not possible, a choked event is formed, which 
could in principle be detected from neutrino observations 
(Meszaros and Waxman, 2001). A distinction is made be- 
tween prompt and delayed collapse to a black hole, with 
black-hole production occurring in the latter case due to de- 
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layed fallback onto a young neutron star (MacFadyen et al., 
2001). 

The collapsar model must contend with the difficulty of 
ejecting baryon-clean material through an overlying shell of 
material. This is accomplished through an active central en- 
gine that persists for at least as long as the prompt phase of 
the GRB. Formation of relativistic jets of baryonic-clean ma- 
terial reaching Tq <~ 10 2 -10 3 represents a major difficulty in 
these models (Tan et al., 2001). 

4.3 Supranova Model 

A central motivation of the supranova model of Vietri and 
Stella (1998) is to identify a site that is originally free of 
baryon contamination. This occurs through a two-step col- 
lapse to a black hole, where a "supramassive" neutron star 
(i.e., with mass exceeding the Chandrasekhar limit) is formed 
in the first-step through a supernova explosion. The neutron 
star is initially stabilized against collapse by rotation. The 
loss of angular momentum support through magnetic dipole 
and gravitational radiation leads to collapse to a black hole 
after some months to years. A two-step collapse process 
means that the neutron star is surrounded by a supernova 
shell of enriched material which can explain rebrightening 
events, as seen in GRB 970508 (Vietri et al., 1999). Alter- 
nately, the neutron star could be driven to collapse by accret- 
ing matter in a binary system (Vietri and Stella, 1999). 

The period of activity of a highly magnetized neutron star 
preceding its collapse to a black hole can produce a pulsar 
wind bubble consisting of low density, highly magnetized 
pairs (Konigl and Granot, 2001), in accord with afterglow 
model fits (Panaitescu and Kumar, 2001). The earlier super- 
nova could yield ~ 1M Q of Fe in the surrounding vicinity. 
The discovery of variable Fe absorption in GRB 990705 dur- 
ing the prompt emission phase (Amati et al., 2000) and X-ray 
emission features in the afterglow spectra of GRB 991216 
(Piro et al., 2000) provide a test of this model. Variable X- 
ray absorption is essentially undetectable from a GRB sur- 
rounded by a uniform CBM, even if highly enriched (Bottcher 
et al., 1999). Extremely clumped ejecta (Bottcher et al., 2002) 
or blueshifted resonance scattering in a high-velocity outflow 
(Lazzati et al., 2001) can account for the variable absorption. 

Fig. 9 illustrates the allowed parameter space in the supra- 
nova model, where / refers to the volume filling factor of 
the ejecta. By contrast, collapsar progenitors must produce 
^ O.IM0 of iron within <~ 1 pc, but massive stars that have 
metal rich winds lose too much mass to form a collapsar 
(Fryer, 2001). A He-merger model, where a binary compact 
object merges with the He core of a companion, requires an 
unreasonably short merger timescale (Bottcher et al., 2002). 

4.4 Other Models 

These three classes of models hardly exhaust the possibili- 
ties. Some alternative suggestions include a strong-field mil- 
lisecond pulsar origin, where the collapse of a white dwarf 
to a neutron star produces an <~ 10 15 Gauss magnetar with a 



Supranova 



10* 



10* 



Quick 
He-merger 



He-merger 



filfiiips 

ion 



Standard 
-mc 

I ; 



krf\ 




10 J 



Fig. 9. Parameter space giving the amount of supernova ejecta mass 
M c i that must be concentrated in clumps versus the time delay At 
between the primary supernova explosion and the subsequent GRB 
to explain the variable absorption feature in GRB 990705 (Bottcher 
et al., 2002). Solid lines give the probability P l c ° a that a clump will 
be along the line of sight, which depends on the iron enhancement 
XFe with respect to Solar system values. The vertical dotted lines 
refer to GRB blast wave Lorentz factors r = 1000r3 that give 
results consistent with the observations. 



millisecond period (Usov, 1992, 1994). Others include the 
"cannonball" model (Dar and De Rujula, 2000) involving 
highly relativistic and highly beamed jets of materials from 
SNe. Another involves the transition of a neutron star in a 
low-mass binary X-ray system to a strange star (Cheng and 
Dai, 1996). In others, a pair electromagnetic pulse (Ruffini 
et al., 2001) or the Blandford-Znajek effect (Lee et al., 2000) 
extract the energy during the black-hole formation event. 

4.5 Best-Bet Model 

In light of the Fe line observations and the constant energy 
reservoir result of Frail et al. (2001), the supranova model 
seems to offer certain advantages over a collapsar model to 
explain the origin of GRBs. For example, 

1 . The narrow distribution of energy releases follows from 
a progenitor model involving neutron stars that collapse 
to black holes, since the neutron-star progenitor masses 
must be very similar. By contrast, a large diversity of 
core masses is possible in the collapsar model, and the 
varying mantle structure would strongly modulate the 
energy radiated as 7 rays. 

2. The immediate environment is much cleaner in the vicin- 
ity of a collapsing neutron star than in a massive-star 
collapsar model. Thus the baryon-contamination prob- 
lem is more easily solved in the supranova model. 

3. The production of large amounts of iron in the GRB 
environment is possible in a two-step collapse process, 
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with the delay between the first- and second-step event 
permitting the nickel and cobalt to decay to iron prior to 
the GRB. 

4. The lack of evidence for winds from massive stars and 
the low surrounding densities is contrary to expectations 
in a model where GRBs originate in the direct collapse 
of massive stars to black holes (Panaitescu and Kumar, 
2001). Moreover, the low densities inferred for the sur- 
roundings are much less than expected for typical ISM 
environment, and could be expected in pulsar wind bub- 
bles (Konigl and Granot, 2001). 

5. The typical E p k energies observed in BATSE can be 
explained in a prompt external shock model, which is 
compatible with the prompt collapse of a neutron star to 
a black hole in the supranova model. The narrow distri- 
bution of Epk energies seems to require fine tuning in 
a model with internal shocks in a relativistic wind, and 
no argument from first principles has been advanced to 
explain the typical durations of the long-duration GRBs 
in a collapsar wind model. 

The evidence is not yet compelling to rule out either model; 
for example, iron line observations can be explained in the 
context of the collapsar model through a persistent source of 
emission involving either a magnetar or the delayed accretion 
onto a black hole (Rees and Meszaros, 2000; Meszaros and 
Rees, 2001). By contrast, the delay between neutron star for- 
mation and collapse to a black hole in the supranova model 
seems to require ad hoc assumptions about the spin rates and 
magnetic fields of the neutron stars formed in SNe explo- 
sions. Various anisotropic and clumpy ejecta are required. If 
GRB history is any guide, the best bet is on a model yet to be 
proposed. 

5 Gamma Ray Burst Cosmology 

We mention some of the main lines of research in GRB cos- 
mology, making no claims for completeness. The rate of 
lensing events on GRBs can set limits to the cosmological 
abundance of dark matter in the form of compact objects, 
depending on the specific cosmology (Marani et al., 1999). 
High-redshift GRBs provide probes of the early universe that 
can be used to examine the redshift dependence of metallic- 
ity, the onset of the first generation of stars, the evolution 
of Lya forest and metal absorption line systems, and the 
epoch of reionization (Lamb and Reichart, 2000; Ciardi and 
Loeb, 2000). An interesting feature of the afterglow is that its 
brightness is not strongly dependent on redshift. The metal- 
licity also may depend on the location of the GRB source in 
the host galaxy (Ramirez-Ruiz et al., 2002). X-ray obscura- 
tion of high-redshift GRBs and the fraction of optically dark 
bursts will reveal whether GRBs are hosted in luminous in- 
frared galaxies similar to the high-redshift SCUBA galaxies 
(Ramirez-Ruiz et al., 2002a). 

Totani (1997, 1999) and Wijers et al. (1998) showed that 
if GRB sources have a massive star origin, then the redshift 



distribution of GRBs yields the star formation rate (SFR) his- 
tory of the universe. Gamma rays are less attenuated during 
transport than optical radiation, so that a complete sample of 
GRB redshifts will arguably give a better measure of the SFR 
than measurements inferred from optical and UV surveys of 
galaxies, which are subject to strong extinction corrections. 

Before a large number of GRB redshifts were known, it 
was first necessary to determine whether the BATSE count 
rate distribution was in accord with the assumption that the 
rate of GRBs followed the SFR history, as inferred, for ex- 
ample, from HST observations (Madau et al., 1998). This 
depends sensitively on the assumed luminosity function of 
GRBs. Krumholz et al. (1998) demonstrated, however, that a 
variety of luminosity and redshift distributions of cosmologi- 
cal GRBs are consistent with the observed peak flux distribu- 
tion, so that a GRB rate density distribution that followed the 
SFR of the universe was not unique. The degeneracy of the 
unknown GRB luminosity function can be broken by assum- 
ing a specific model for GRB emission, as done by Bottcher 
and Dermer (2000) for an external shock model. This leads 
to a specific prediction for the redshift distribution of GRBs 
(Fig. 2, lower panel) based upon the assumed SFR function, 
and the GRB energy and Lorentz factor distributions. Com- 
parison with GRB redshift data can then be used to refine the 
model parameters and test the basic assumption that the rate 
density of GRBs follows the SFR of the universe. 

GRBs can have important effects on the galaxy environ- 
ment, and were proposed as the energy sources of HI shells 
and stellar arcs (Efremov et al., 1998; Loeb and Perna, 1998), 
though the evidence for beaming makes these possibilities 
less likely. GRBs could melt dust grains by GRB UV radi- 
ation to produce flash-heated chondrules in the early Solar 
system (McBreen and Hanlon, 1999). In one sense, the most 
cosmologically interesting effect of GRBs is their influence 
on biological activity (Scalo and Wheeler, 2002) and life ex- 
tinction (Dar and De Rujula, 2001), which depends not only 
on GRB photon production but on particles accelerated by 
GRBs, namely cosmic rays. 

6 Cosmic Ray Production by GRBs 

As summarized in §4, several lines of evidence indicate that 
GRBs are closely related to a subset of SNe that drives a 
relativistic outflow in addition to the nonrelativistic ejecta 
expelled during the collapse of the massive core to a neu- 
tron star. The relativistic ejecta decelerate to nonrelativis- 
tic speeds at the Sedov radius by sweeping up matter from 
the external medium. Particle acceleration occurs at these 
shocks, just as in the nonrelativistic shocks from "normal" 
Type la and Type II supernova remnants, though with the ad- 
dition of a relativistic phase of deceleration. The first-order 
shock Fermi mechanism is generally recognized as the mech- 
anism that accelerates GeV/nuc cosmic rays in the converg- 
ing flows formed by the forward shock. In addition, second- 
order Fermi acceleration of particles through gyroresonant 
particle-wave interactions with magnetic turbulence gener- 
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ated in the shocked fluid can also accelerate particles to high 
energies. 

We consider the simplest case of particle acceleration via 
the first- and second-order Fermi processes at the external 
shock formed by ejecta decelerating in a uniform surround- 
ing medium. These considerations will establish whether 
GRBs can accelerate particle to the ankle of the cosmic ray 
spectrum at 3 x 10 18 -10 19 eV and to ultra-high (> 10 19 eV) 
energies. We show that ultra-high energy cosmic rays (UHE- 
CRs) can in principle be accelerated by GRB blast waves. 
Additional effects involving internal shocks, acceleration at 
the reverse shock, or a highly magnetized CBM can improve 
the acceleration efficiency. 

6. 1 Constraints on Acceleration 

A basic requirement for particle acceleration to some maxi- 
mum energy E max in the observer frame is that the Larmor 
radius r' L be smaller than the comoving size scale r' of the 
system. The comoving Larmor radius r' L = mc 2 j/qB for 
particles with mass m = Am p and charge q = ze. The 
requirement that r' L < r' implies that 
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eV , (23) 



where r' is related to the measured variability timescale t var 
through the relation t var = (l+z)t' var /5 = ( 1 +z) r'/c5, and 
5 w r. Recalling from equation (9) that B = 116^/eBno^3oa 
G, we see that UHECR acceleration can be achieved for highly 
relativistic GRB outflows, especially if large Z nuclei such as 
iron are accelerated. 

In Fermi processes, particles can increase their energy by 
<~ mc 2 7 on a Larmor timescale = = 2irmc^ / qB 

(Rachen and Mesazros, 1998). Hence the maximum accel- 
eration rate is j acc = f acc (qB/2wmc), where f acc < 1. A 
basic constraint on particle acceleration is the radiation re- 
action limit due to synchrotron losses. The synchrotron loss 
rate 
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Equating this with A / acc gives the radiation-reaction limit on 
particle acceleration through Fermi processes 
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This gives E maXtSyn w 6 x 10 15 r 300 \/facc/B(G) eV for 
electrons, and a maximum synchrotron photon energy hv max 
= 20/ acc r MeV, independent of B, as is well-known (de 
Jageretal., 1996). For ions, the synchrotron radiation-reaction 
limiting energy is 
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The maximum ion synchrotron photon energy is hv max = 
40(j4/ acc /Z 2 )rGeV. Other processes, such as Compton losses 



for leptons and photomeson and photopair losses for ions, 
can limit acceleration further. Provided that these losses are 
not severe, the synchrotron radiation reaction is not a severe 
limitation for UHECR particle acceleration. 

Another kinematic limitation on first-order Fermi acceler- 
ation arises from the size scale of the system. Because par- 
ticles must be scattered upstream of the shock, the upstream 
Larmor radius tl must be smaller than the characteristic size 
of the system during the period when the shock is strong. 
The Sedov radius, given by equation (5), is the radius where 
the blast wave has swept up an amount of matter equal to 
the ejecta mass. Particle acceleration is strongly suppressed 
after the explosion enters the Sedov regime. The condition 
r*L < is implies that 
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(27) 



where the B^q and n are the mean magnetic field, in /j,G, 
and density of the interstellar medium. 

The time required to cycle across the supernova shock pro- 
vides a further limitation on particle acceleration by shock 
Fermi processes. The first-order Fermi acceleration rate 'Jfi 
w P^/tcyc, where is the shock speed, and the cycle time 
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Here k is the diffusion coefficient, u is the speed of the flow 
in the rest frame of the shock, and subscripts — and + re- 
fer to quantities upstream and downstream of the shock, re- 
spectively (e.g., Drury, 1983; Kirk, 1994). The expression 
on the right-hand-side assumes that the time to scatter is de- 
termined by quantities upstream of the shock, and that the 
diffusion coefficient is given by the Bohm diffusion limit 
nrL,-. Differences in cycle times from quasi-parallel and 
quasi-perpendicular shocks arise from shock-drift in the lat- 
ter case. The maximum energy-gain rates per unit distance 
x for nonrelativistic quasi-parallel and quasi-perpendicular 
shocks, respectively, are given by 
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(29) 



The energy increases by only a factor = 2 in a single cy- 
cle following the first cycle in relativistic shock acceleration 
(Gallant and Achterberg, 1999). 

By integrating the energy-gain rate over the equation of 
blast wave evolution given by equation (6), the classic re- 
sult of Lagage and Cesarsky (1983) is generalized to give the 
maximum particle energy 



E„ 



~in 16 7R r 2 / z ( !^©£o 
no 



) 1 / 3 eV 



(30) 



that can be obtained by first-order Fermi acceleration at a 
forward external shock. This shows that neither nonrelativis- 
tic or relativistic first-order shock-Fermi processes are capa- 
ble of accelerating particles to the ankle of the cosmic ray 
spectrum. In fact, nonrelativistic (0o^o *C 1) shocks have 
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Fig. 10. Maximum proton energies from first- and second-order 
Fermi acceleration for explosions with initial Lorentz factors To 
and speed (Hoc. 

difficulties to accelerate particles past the knee of the spec- 
trum (Lagage and Cesarsky, 1983), although winds from the 
pre-supernova star can enhance the local magnetic field and 
reduce the surrounding medium density and so permit accel- 
eration to much higher energies (Volk and Biermann, 1988). 

Second-order Fermi acceleration from gyroresonant inter- 
actions between particles and waves in relativistic shocks can 
accelerate particles to considerably higher energies than first- 
order Fermi acceleration when the wave speeds approach c. 
For a wave spectrum composed of Alfven waves with speed 
v\, the acceleration rate from a power law spectrum of waves 
w(k) oc k~ v with wavenumber k > k min is given by 

iFii « ( — - 1) (^(rikminiy- 1 ■ (3D 
c rl 

Here r° = mc 2 /qB is the nonrelativistic Larmor radius, and 
v is the index of the turbulence spectrum (v — 5/3 corre- 
sponds to a Kolmogorov spectrum). Magnetic field turbu- 
lence is produced in the relativistic shocked fluid when parti- 
cles and dust are captured and isotropized (Pohl and Schlick- 
eiser, 2000; Schlickeiser and Dermer, 2000). After integrat- 
ing over the energy gain rate during the duration of a decel- 
erating blast wave, one obtains a maximum particle energy 

E maxJ i w 8 x lO 2O Zi^4 /2 nJ /6 /A/3o(™ ro) 1/3 eV (32) 

(Dermer, 2001), where K v ^ [Q£e B MI f&\ 1/{2 ~ v) forv = 
3/2,5/3, and f&x/T 2 ^ x/\2T 2 gives the observer frame 
width of the blast wave from the shock-jump conditions in a 
uniform CBM (Vietri, 1998, 1995; Panaitescu et al., 1998; 
Dermer and Humi, 2001). The strong dependence of the 
maximum energy on f3 makes second-order processes unim- 
portant in nonrelativistic flows when (3 <C 0.1. 

Fig. 10 compares the maximum proton energies that can 
be obtained by first-and second-order Fermi acceleration at 
an external shock formed by an explosion with total energy 



m = 1, corresponding to 1.8 x 10 5 ergs. The surround- 
ing external medium contains a 1 /iG field and a density 
of 1 cm~ 3 . A fully turbulent magnetic field with £ = 1, 
cb = 0.2, and v = 3/2 is assumed in the second-order pro- 
cess. First-order Fermi processes cannot accelerate particles 
to ultra-high energies for particles accelerated from low en- 
ergies for typical ISM conditions. Relativistic shocks can, 
however, give a pre-existing relativistic particle population a 
single-cycle boost by a factor T 2 (Gallant and Achterberg, 
1999; Vietri, 1995). 

6.2 UHECR by GRBs 

The Larmor radius of a particle with energy 10 20 i?2o eV is 
w 100 E 20 /(ZB tiG ) kpc. Unless UHECRs are heavy nu- 
clei, it is unlikely that they originate from our Galaxy. The 
proposal that UHECRs are accelerated by extragalactic GRB 
sources was first advanced by Vietri (1995), Waxman (1995), 
and Milgrom and Usov (1995). The first two of these authors 
noted that the energy density of UHECRs was comparable to 
the energy density that would be produced by GRB sources 
within the GZK radius, assuming that the 7-ray and UHECR 
power from GRBs are roughly equal. Although the details 
of this estimate are incorrect inasmuch as they were made 
before the Beppo-SAX results on the detailed redshift dis- 
tribution of GRBs, the overall energy-density estimate still 
gives strong circumstantial support for this hypothesis. 

The argument proceeds by noting that the energy density 
of UHECRs observed near Earth is simply 

> LcRBtesc ,,,,,, 

UUH = C —r (33) 

Vprod 

where Lgrb is the power of GRBs throughout the produc- 
tion volume V pro d of the universe. UHECRs are produced 
with an efficiency £ compared with the 7-ray power and "es- 
cape" from the universe primarily through photohadronic pro- 
cesses with an effective escape time t esc . 

The mean 7-ray fluence of BATSE GRBs is F 1 w 3 x 
10~ 6 ergs cm~ 3 and their rate over the full sky is Ngrb ~ 
2/day. If most GRBs are at redshift (z) <~ 1, then their 
mean distance is (d) s=s 2 x 10 28 cm, so that the average 
isotropic energy release of a typical GRB source is (E 7 ) 
47r(«i) 2 F 7 /(l |z)=8x 10 51 ergs, implying a mean GRB 
power into the universe of Lgrb ~ 2 x 10 47 ergs s _1 . 
(This estimate is independent of the beaming fraction, be- 
cause a smaller beaming fraction implies a proportionately 
larger number of sources.) UHECR protons lose energy due 
to photomeson processes with CMB photons in the reaction 
p + 7^p + 7r°,n + 7r + . The effective distance for 10 20 eV 
protons to lose 50% of their energy is 140 Mpc (Stanev et al., 
2000), so that t esc = 1.5 x 10 16 . This implies that the energy 
density observed locally is u UH ~ 10~ 22 ri ergs cm~ 3 . By 
comparing with Fig. 1 1, we see that this estimate is about an 
order of magnitude below the energy density of super-GZK 
particles with energies exceeding 10 20 eV. 

Given the limitations of this crude estimate, the coinci- 
dence between the GRB power and the UHECR energy den- 
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Fig. 11. Cosmic ray energy density observed near Earth. The 
two solid curves show extrapolations to model fits of the observed 
cosmic-ray proton (Simpson, 1983) and all-particle (Fowler et al., 
2000) spectrum, and the data points are UHECR observations from 
AGASA (Takeda et al., 1998). 



sity suggests that this hypothesis may be correct if GRBs can 
accelerate UHECRs with high efficiency (t] w 1). An im- 
proved estimate must consider several other aspects, includ- 
ing 

1 . The efficiency of GRB sources to produce 7-ray emis- 
sion within the BATSE band. If only 10% of the to- 
tal energy of GRBs are radiated at soft 7-ray energies, 
then these sources will be much more energetic on aver- 
age, which will improve the comparison - provided that 
UHECR acceleration remains very efficient. 

2. The actual redshift distribution of GRB sources. If the 
rate density of GRB sources follows the SFR history of 
the universe (Madau et al., 1998), then the local comov- 
ing rate density will be w 6-10 times less than the rate 
density at z w 1. This effect, when factored into our 
improved knowledge of the GRB redshift distribution 
showing that the faintest GRBs are at z > 1, makes the 
hypothesis that UHECRs are accelerated by GRBs more 
difficult to sustain (Stecker, 2000). 

3. The existence of the dirty and clean fireball classes of 
GRBs. These sources classes would not trigger BATSE, 
and the dirty fireball class could be detected as the X-ray 
rich GRBs (Heise et al., 2001). This correction could 
increase the average power into the universe by a factor 
- 3 - 5 (Dermer, 2000). 

4. The uncertainty in the energy density of the super-GZK 
cosmic rays. The UHECR energy density shown in Fig. 
1 1 is from the AGASA measurements. Monocular HiRes 
observations show a lower flux than AGASA above w 
10 20 eV by a factor of w 5 (Sommers, 2002), which 
improves the coincidence between the UHECR energy 
density and the GRB power. 



A detailed estimate of GRB power is given by Dermer (2000) 
in the context of an external shock model for GRBs per- 
formed by Bottcher and Dermer (2000). The conclusion of 
this study is that GRBs remain a viable source origin for 
UHECRs. A testable prediction of the hypothesis is that star- 
forming galaxies which host GRB activity will be surrounded 
by neutron-decay halos (Dermer, 2000), as discussed in §7. 

6.3 Rate and Power of GRBs in the Milky Way 

The preceding estimates were made for extragalactic sources 
of cosmic rays. GRBs will also take place in galaxies such 
as our own. After a relativistic GRB blast wave decelerates 
to nonrelativistic speeds, it evolves in a manner similar to 
nonrelativistic supernova remnant shocks. GRBs will there- 
fore power cosmic rays with energies below the knee of the 
cosmic ray spectrum. These cosmic rays will add to cosmic 
rays accelerated in the remnants of Type I and Type II SNe. 
Here we estimate the rate and power of GRBs throughout the 
Milky Way in order to determine the contribution of GRBs 
to cosmic ray production. 

We assume the constant energy reservoir result of Frail 
et al. (2001) to estimate the rate of GRBs into L* galax- 
ies such as the Milky Way. The local density of L* galax- 
ies can be derived from the Schechter luminosity function, 
and is w 1/(200-500 Mpc 3 ) (Wijers et al., 1998; Scalo and 
Wheeler, 2002; Dermer, 2000). The BATSE observations 
imply, as already noted, ~ 2 GRBs/day over the full sky. 
Due to beaming, this rate is increased by a factor of 5OO/500, 
where / 50 o ~ 1. Given that the volume of the universe is 
~ 47r(4000 Mpc) 3 /3, this implies a rate per L* galaxy of 
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(34) 



The factor SFR corrects for the star-formation activity at 
the present epoch. The SFR rate per comoving volume is 
about 1/6 as active now as it was at z = 1, and the factor 
Kpx accounts for dirty and clean fireball transients that are 
not detected as GRBs. Thus a GRB occurs about once every 
5000 years throughout the Milky Way. 

Table 1 gives the outflow speeds and rates of different 
types of SNe and GRBs in the Galaxy. The data for the 
outflow speeds of Types la, II, and Ib/c come from Lozin- 
skaya (1992) and Weiler et al. (2000). The mean values of 
the initial dimensionless momentum (@oTq) of the outflows 
are also given. The supernova rate data are obtained by mul- 
tiplying the results of Cappellaro et al. (1999) for galaxies 
of Type Sbc-Sd (renormalized by Panagia (2000) to an H- 
band luminosity for h 7a = 1) by a factor of 2 to account 
for the luminosity of the Milky Way. One sees that the GRB 
rate is about 10% as frequent as Type Ib/c SNe, and about 
1% as frequent as Type II SNe. The rates of dirty and clean 
fireballs are unknown, but if the X-ray rich 7-ray bursts com- 
prise the tail of the dirty fireball population, then they could 
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Table 1. Ejecta Speeds and Rates of Supernovae and GRBs in the 
Milky Way 



Explosion 


Outflow 


(polo) 


Rate 


Type 


Speed (km s _1 ) 




(century ~ 4 ) 


SN la 


< 2 x 10 4 


0.03 


0.42 


SNII 


- 10 3 -2 x 10 4 


0.01 


1.7 


SN Ib/c 


~ 1.5 x 10 3 -2 x 10 5 


0.2 


0.28 


Dirty 


3 x 10 s 


30 


? 


GRB 


3 x 10 s 


300 


~ 0.02 


Clean 


3 x 10 s 


3000 


? 



be even more frequent than the GRB population. Statistical 
fits (Bottcher and Dermer, 2000) to the BATSE data suggest 
that the rate of clean fireballs is smaller than the GRB rate. 

By weighting the GRB rate, equation (34), by the mean en- 
ergy of 3 x 10 51 / (?7 7 /0.2) ergs, we see that the time-averaged 
power of GRBs throughout the Milky Way or other L* galax- 
ies is 



2 x 10 40 



f: 



500 



(77 7 /0.2) v 1/6 



SFR Kft, 



•(35) 



Because each GRB has a total energy release of a few x 10 51 
ergs, comparable to SNe, but occur ~ 100 times less fre- 
quently, the available power from the sources of GRBs is 
only a few per cent of the power from Type la and Type II 
SNe. This would suggest that GRB sources can only make a 
minor contribution to cosmic ray production below the ankle, 
but the efficiency for accelerating cosmic rays in the nonrel- 
ativistic SN outflows could be considerably less than in the 
relativistic outflows of GRBs. 

In this regard, note that Cas A has recently been detected 
with HEGRA at TeV energies at the level of w 8 x 10~ 13 ergs 
cm~ 2 s _1 (Aharonian et al., 2001). Spectral considerations 
suggest that the emission is hadronic. The expected level of 
100 MeV- TeV 7-ray production can be estimated to be 

F 7 w (47rrf 2 ) -1 x 10 51 ergs rj p x c<7 pp n « 



7 x 10 



d 



3.4 kpc 



~ 2 ergs cm 2 s 



(36) 



where the strong interaction cross section a vp k, 30 mb. This 
expression may underestimate the spectral energy flux mea- 
sured at a fixed frequency range and the 7-ray production 
efficiency, but unless the density is unusually tenuous in the 
vicinity of Cas A and other SNRs examined with the Whip- 
ple telescope (Buckley et al., 1998), the TeV observations 
suggest that most SNRs do not accelerate cosmic rays with 
high efficiency. 

The hypothesis that cosmic rays are accelerated by SNRs 
is favored in view of the power provided by SNe (Gaisser, 
1990). With an energy release of ~ 10 51 ergs per SN and a 
rate of w 1 every 30 years, the power into the Galaxy from 



SNe is at the level of sa 10 42 ergs s _1 , well in excess of the 
« 5x 10 40 ergs s _1 required to power GeV/nuc cosmic rays. 
Moreover, first-order Fermi acceleration by a strong shock 
produces a —2 spectrum, similar to the injection spectrum 
required to explain the measured cosmic ray spectrum prior 
to steepening due to propagation effects. 

In spite of this apparent success, observational difficulties 
for the hypothesis that cosmic rays are accelerated by SNe 
bear repeating (Dermer, 2000) 

- The unidentified EGRET 7-ray sources have not yet been 
firmly associated with SNRs and do not display distinct 
7T° features (Esposito et al., 1996); 

- The spectrum of the diffuse galactic 7-ray background 
is harder than is expected if the locally observed cosmic 
ray proton and ion spectrum is similar throughout the 
Galaxy (Hunter et al., 1997); 

- The acceleration of particles above the knee of the cos- 
mic ray spectrum is difficult to explain in the simplest 
theory (Lagage and Cesarsky, 1983) of cosmic ray ac- 
celeration by SNRs; and 

- As we have seen, TeV radiation is not detected at the 
expected levels. 

In view of these difficulties, it seems possible that cos- 
mic rays are predominantly accelerated by the subset of SNe 
which eject relativistic outflows. Bright enhancements of 
hadronic emission from those SNe which host GRB events 
are implied. From the rate estimates shown in Table 1, we 
see that about 1 in 20 to 1 in 100 SNRs would exhibit this 
enhanced emission. The better imaging and sensitivity of the 
GLAST telescope and the next generation of ground-based air 
Cherenkov telescopes will be able to test this hypothesis. 



7 High Energy Neutrons and Neutrinos from GRBs 

Observations of neutral particles and their secondary radia- 
tions will test the hypothesis that cosmic rays are acceler- 
ated by GRB blast waves. Neutron-decay halos are formed 
around galaxies that harbor GRB activity (Dermer, 2000). 
Equation (35) shows that the sources of GRBs inject a mean 
power of 10 40±1 ergs s~ 4 into a typical L* galaxy. In order 
to explain the locally observed UHECR energy density, each 
GRB source must accelerate ions with high efficiency. Neu- 
trons will be directly produced through the process p + 7 — > 
n + 7T+ in the GRB blast wave and will escape the accelera- 
tion site. The neutrons, unbound by the magnetic field in the 
blast wave, leave the acceleration site with Lorentz factors 
7„ = 10 10 7io, where 0.1 £ 710 £ 100. The neutrons decay 
on a timescale j n t n = 3 x 10 5 7io yrs, where the neutron 
/3-decay lifetime t n = 900 s. 

The neutrons travel a characteristic distance A„ = 9O710 
kpc before they decay into highly relativistic electrons and 
protons. The electrons initiate an electromagnetic cascade 
that forms a diffuse w 10 TeV Compton 7-ray component 
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and a nonthermal synchrotron component extending to soft 
7-ray energies. The halo is unfortunately very weak. Three 
orders of magnitude reduction in the injection luminosity re- 
sults from The /3-decay to an electron results in three orders 
of magnitude reduction of power, and a further reduction re- 
sults from the efficiency of photomeson production in the 
GRB blast wave. In an external shock model, the efficiency 
can be lower than 1 %, implying powers of neutron decay ra- 
diation halos surrounding L* galaxies at <~ 10 35 ergs s _1 . 

Will water and ice Cherenkov neutrino telescopes detect 
GRBs in the light of high-energy photohadronic neutrinos? 
The comoving photon energy density is 



eV(e') 



4 

2cr> 2 



^2 , and e : 



2Te' 
1 + z 



(37) 



for a blast-wave geometry. The proper-frame density increases 
by 4T upon crossing a relativistic shock. The invariance of 
the total number of particles implies that/y = x 2 (Ail)A' sh n' 
= (Afl)n f* dxx 2 , so that 



r' = A' = - 

b sh 12Yx 2 



jfp , uniform CBM 



J-i'D .r-Tg „ i 

4r ' x ^ 



(38) 



The comoving photon density is n'(e') = e'u'(e')/m e c 2 e' 2 , 
and the synchrotron emission is assumed to be isotropic in 
the comoving frame, so that n'(e', //) = ^n'(e'). 

The time scale for significant energy loss by photohadronic 
reactions is 

c r°° f 1 

t' p ~^ ~ - jf de>J i df,'(l- t ,')n'(e , , l i')a PJ ^(e") , (39) 

where a PJ ^ w (e") = a a S[/j,' - (1 - e A /^' p e')], a w 2 x 
1(T 28 cm 2 , e A « 640, e" = ~/' p e'(l - /u')> and we now use 
primes to refer particle Lorentz factors to their proper frame. 
The comoving time available to undergo hadronic reactions 
is t ava = r' b /c = 2Yt/ (1 + z). The quantity 



a d 2 L e A 



20%m e c 3 r 2 r' b 



de' e'" 3 /, 



(40) 



represents an efficiency for converting accelerated particle 
energy into neutrinos. The threshold condition ~f' p e' = ca re- 
quires acceleration of protons with observer-frame energies 
7 P = r 2 e A /[(l + z)e pk ] oc T~ 2 £ i 3 / 4 from scattering with 
photons at the vF v peak photon energy e p k- As can be seen, 
demands on particle acceleration are more easily satisfied 
during the prompt phase of the GRB. Neutrinos are formed 
in the observer's frame with energies of w 0.05to p c 2 7p . The 
peak photon energy of the vF v spectrum is 



2T 



€pk 



(1 + *) 



500 

l + z 



^oo- (41) 
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The brightest neutrino production episode occurs t ~ td for 
a fast cooling scenario, as can be shown from equation (40). 



Quantities evaluated at t d wear hats. Equations (13) and (15) 
or rule (18) (Dermer et al., 1999a) imply 
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(42) 



so / e = fe pk (e/i p k)°' 1 ' and a„ = 1/2 for e < e pk , a v = 
(2 — p)/2 for e > e p fe in the fast cooling limit. From eq. (40), 
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(44) 



letting r£ = f A x d /T. 

The low efficiencies of neutrino production makes detec- 
tion of GRBs with a 1 km 3 neutrino telescope extremely un- 
likely within the context of a uniform external shock scenario 
(Dermer, 2000). Low Lorentz-factor dirty fireballs have in- 
creased efficiencies f]f c by as much as three orders of magni- 
tude, but require acceleration to ^ pk ~ 10 10 /r 2 fc 2 e 2 v / e7i7jo, 
which follows from the relation y pk = r§eA/7 P fe(l+2 ; )epfe = 
1. The fast-cooling condition -f c < jmin implies e e es ^ 
4x 10~ 5 /[El^nl^kpil+z)^}, which additionally shrinks 
available parameter space. 

Result (44) gives the photomeson production efficiency at 
the deceleration time scale for a proton with energy nip^jph 
~ 10 14 eV/[fc 2 e 2 ^esn r 2 00 ] eV in a fast cooling blast 
wave. Secondary neutrinos with energy <~ 15% xm p c 2 -f p k 
would be received from such GRBs. Deep water (Baikal, 
Nestor, Nemo, Antares) and ice (AMANDA, IceCube) neu- 
trino detectors have detection probability P <~ 10~ 4 P_4, 
with P_4 ~ 1 at PeV neutrino energies. The number of de- 
tected neutrinos N v can be approximated by the product of 
the telescope area 10 10 A w cm 2 , the neutrino flux = 10 52 -E 52 
/[47rd 2 : (mpC 2 7pfe)] cm~ 2 s _1 , the multiplicity 3/2 formuon 
neutrinos, and a bandwidth correction factor 0.1£_i. The 
threshold energy for photomeson production by protons in- 
teracting with photons at the peak of the vF v synchrotron 
spectrum at t = td is given by m p c 2 -f p k — w p c 2 rgeA/[(l + 
z)i P k] = 170 eigs/[k 2 e 2 e ^/eWn^rl m }. 

Multiplying these factors together gives an estimate of the 
number of neutrinos that would be detected from a GRB, as- 
suming that a comparable amount of energy is in the non- 
thermal proton spectrum as in the explosion: 



JV„~4x 10" 



AioP_ 4 £^ 3 fcpe e no /3 r3^Co.i 
d 2 s (2-a„, u )(l + z)f A 



(45) 



Equation (45) is consistent with the numerical calculations 
of Dermer (2000), but much less accurate by neglecting the 
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energy dependence of P, and neutrino production in the af- 
terglow phase. This result shows that there is no prospect 
for detecting £ 1 neutrinos from FRED-type GRBs, and we 
verify the computations (Dermer, 2000) of the weak neutrino 
flux produced by GRB blast waves decelerating in a uniform 
CBM. 

Equation (45) says that GRBs are much brighter v sources 
when r 30 o 3> 1, n 3> 1, and /a <C 1. Clean fireballs 
are however rare (Bottcher and Dermer, 2000). Density in- 
homogeneities from the SN envelope in the supranova model 
could enhance v production in an external shock model, not 
only through the density effect indicated in equation (45), 
but through photo-hadronic processes with external photons 
(Atoyan and Dermer, 2001) and back-scattered photons. Eq. 
(38) shows that thin shells are formed in a clumpy medium, 
so that /a could be <C 1/12. 1 Because spiky and highly vari- 
able GRBs are due to large density contrasts in the external 
shock model on a size scale <C x/T (Dermer and Mitman, 
1999), it seems possible that the spikiest and most variable 
GRBs could produce detectable high-energy v emission in 
an external shock model. A definitive statement must await 
further studies. 

Internal shock-model advocates (Waxman, 1995; Waxman 
and Bahcall, 1997) determine the blast-wave width from the 
observed variability time scale t var ~ (1 + z)r' b / (cS), so that 
r' b = /ax/T = 2cTt var /{\ + z). Thus / A ~ t var /t d , that 
is, the shell must lie within a distance ct var <C ctd of the 
central source, where the deceleration time is associated 
with any surrounding CBM. Sufficiently small shell widths, 
though larger than the size scale of the original explosion, 
can enhance v production considerably, so that the observ- 
able 7-ray power sets the only available scale on energy re- 
lease. Pair attenuation in these inner regions will make 7-7 
spectral cutoffs in the GLAST energy range. Studies of 7-7 
opacity during the prompt phase of an external shock model 
find transparency up to TeV energies (Dermer et al., 2000a). 

In the uniform external shock model, as we have seen, a 
very low efficiency for neutrino production is a consequence 
of the low radiation energy density in the comoving frame 
(Dermer, 2000). In models involving internal shocks, by con- 
trast, the neutrino production efficiency can be much higher 
due to the larger radiation density (Waxman and Bahcall, 
1997, 2000). In the latter case, GRBs should be detectable 
with km 3 scale neutrino detectors if UHECRs are indeed ac- 
celerated by GRB blast waves. Photon emission from the re- 
verse shock can also enhance neutrino emission (Dai and Lu, 
2001). Neutron decoupling from protons in the expanding 
fireball can also produce a strong flux of 5-10 GeV neutrinos 
(Bahcall and Meszaros, 2000). Secondary production can 
also be effective for high-energy neutrino production (Schus- 
ter et al., 2001). 

We think that no measurable neutrino flux will be detected 
from smooth-profile GRBs. Spiky GRBs could be neutrino- 
bright, that is, potentially detectable in a km 3 -scale experi- 

'This answers a criticism of the external shock model posed by 
R. Sari. 



ment. It will be important to have continuous sky coverage 
with a GRB/hard X-ray survey instrument as IceCube and 
the northern hemisphere water arrays accumulate data. 

8 Advanced Blast Wave Theory 

All aspects of GRB theory cannot be treated in a short re- 
view, but at least some mention should be made of the main 
subjects overlooked. The physics of colliding shells in an 
internal shock/wind scenario was omitted, in keeping with 
Ockam's razor: "pluralitas non est ponenda sine neccesi- 
tate," that is, "entities should not be multiplied unnecessar- 
ily." This hypothesis is not needed in GRB research; for an 
opposite viewpoint, see Piran (1999), and for a balanced as- 
sessment, see Meszaros (2002). (Colliding shell physics is 
of course important in black-hole jet studies; see, e.g., Spada 
et al. (2001).) 

The physical approach summarized here has been devel- 
oped by dozens of authors, and represents a generalization 
of supernova remnant physics to explosions with relativistic 
ejecta. Some important early contributions leading to the rel- 
ativistic synchrotron-shock model include the works of Katz 
(1994) and Tavani (1996). A departure from a nonthermal 
synchrotron origin for the prompt emission is made by Laz- 
zati et al. (2000) and Ghisellini et al. (2000), who argue that 
it is instead due to bulk Compton upscattering in a high pho- 
ton energy-density region in the funnel of an exploding star. 
Part of the motivation for their work was observations of 
very hard X-ray emission with vF v indices a v > 2 during 
the prompt 7-ray luminous phase of 5-10% of GRBs (Crider, 
et al., 1997; Preece et al., 1998). This strongly contradicts the 
optically-thin synchrotron shock model, which predicts that 
only spectra with a v < 4/3 can be radiated in a nonthermal 
synchrotron model. 

Elaborations on the synchrotron-shock theory that may ac- 
count for such hard spectra include photospheric emission 
(Meszaros and Rees, 2000) and radiation reprocessing by a 
medium heated by the 7 rays (Dermer and Bottcher, 2000). 
These latter effects could lead to sites of enhanced annihila- 
tion radiation in the ISM at the location of GRB remnants, 
which INTEGRAL could discover. Madau and Thompson 
(2000) and Thompson and Madau (2000) treated Compton- 
backscatter effects from evolving 7-ray photon fronts, though 
initially approximated as an ^10-100 lt-s shell of photons. 
The Compton back-scattered photons provide targets for suc- 
cessive waves of incident GRB photons through 77 pair- 
production interactions. Beloborodov (2001) and Meszaros 
et al. (2001) consider how pair-loading and preacceleration 
("surfing") of the medium ahead of the 7-ray light front can 
alter the early afterglow and produce a bright optical/UV 
emission component during the early phases of a GRB . These 
early optical flashes might have been mistaken for a reverse 
shock emission. Paczyhski (2001) suggests that these effects 
could produce optical flashes preceding a GRB. 

Finally, we note that the theory of the central engine and 
gravitational radiation from GRBs is not treated here. 
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9 Conclusions 

The Introduction summarizes the scope of this review, and 
will not be repeated here. The themes that I have tried to 
stress, and that would not have been possible to include in a 
review written only a year ago, are the implications of after- 
glow studies and X-ray absorption and line observations on 
source models. In particular, the constant energy reservoir re- 
sult of Frail et al. (2001), derived by modeling breaks in GRB 
afterglow light curves, implies that GRBs are much more nu- 
merous and much less energetic than earlier believed. GRB 
sources are probably related to an unusual type of supernova 
occuring at ~ 1% of the rate of Type II SNe, with energy 
releases several times greater than Type la or Type II SNe. 
Relativistic shock waves can accelerate particles to the ankle 
of the cosmic ray spectrum through second-order Fermi pro- 
cesses, so that cosmic rays between the knee and ankle of the 
cosmic ray spectrum could originate from GRB sources. If 
relativistic shocks accelerate particles with much greater ef- 
ficiency than nonrelativistic shocks, then GRB sources could 
make a significant contribution to the production of cosmic 
rays with energies from GeV/nuc to the knee. 

Cosmic rays may originate from the subset of supernovae 
that eject relativistic outflows. Gamma-ray observations with 
GLAST and ground-based air Cherenkov observatories will 
be fundamentally important to identify the sources of cosmic 
rays by searching for hadronic emission from discrete sites 
in the Galaxy. The Swift Observatory, in conjunction with 
ground-based followup, will establish a large database from 
which GRB redshifts, energies, and beaming fractions will 
be obtained. Neutrino observatories may force GRB source 
models to be revised. We look forward to further advances 
that will make this review obsolete. 
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